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Abstract 

We present a new and novel synthesis of all existing neutrino data regarding the 
disappearance and appearance of Ue and f^. We assume four neutrinos: Vf.-.i-'^i-.i'T-, 
as well as a heavier singlet neutrino Vg of a few eV. The latter may decay into a 
massless Goldstone boson (the singlet Majoron) and a linear combination of the doublet 
antineutrinos. We comment on how this scenario may be verified or falsified in future 
experiments. 



Accepting the totality of present experimental evidence for neutrino oscillations 0, 0, ||, 
it is not unreasonable to entertain the idea that there are four light neutrinos. Since the 
invisible decay of the Z boson tells us that there are only three light doublet neutrinos, i.e. 
z/g, t'^, z/,-, the fourth light neutrino should be a singlet. Usually, is assumed to mix 
with the other neutrinos in a 4 x 4 mass matrix for a phenomenological understanding |Q 
of all the data. However, given that Vg is different from Ve,ii.,Ti it may have some additional 
unusual property, such as decay. In fact, as shown below, this is a natural consequence of 
the spontaneous breakdown of lepton number in the simplest modelf^, and it has some very 
interesting and verifiable predictions in future neutrino experiments. 

If only atmospheric and solar neutrino data are considered, then hierarchical three- 
neutrino oscillations with 



1 
1 

72^ 



v\ = UeCosO 7=(z/^ + z/^) sin6', (1) 

1^2 = J^e sin 9 H — ■j= {u^ + cos 6, (2) 



^^3 = -^(-z/^ + z/^), (3) 

where mi « m2 « m^, would fit the data very well. Here ml ~ 10~^ eV^, (sin^ '20)atm = 1, 
and ~ 10~^ eV^ for the matter-enhanced oscillation solutionis to the solar neutrino deficit 
with (sin^20)soi ~ 10"'^ or near 1, or m"^ ~ 10"^'' eV^ for the vacuum oscillation solution 
with (sin^26')so« ~ 1- 

We now add a fourth neutrino and assume that it mixes a little with z/g and z/^ to 
explain the LSND dataf^. Since the relevant Am^ is now about 1 eV^, it is natural to take 
ml ~ 1 eV^, but this hierarchical solution is disfavored0, because the observed z/^ z/g 
probability^] is contradicted by the z/^ —>■ data of CDHSW|§ together with the z/g Ve 
data of Bugey0. However, there are two ways that this conclusion may be evaded. (1) Let 
777,4 ~ 25 eV^, then the constraint due to the CDHSW experiment is not a factor, but now 



there are three other accelerator v^. experiments: BNL-E734[0, BNL-E776|TT|, and 



CCFR|T^, which have bounds close to but allowed by the LSND 99% likelihood contour. 



This is a marginal hierarchical four-neutrino oscillation solution to all the data. (2) If z/4 
decays, then the parameter space for an acceptable solution should open up. For example, 
in the CDHSW experiment, two detectors at different distances compare their respective Vy, 
fluxes and the ratio is taken. If the component of decays away already before reaching 
the first detector, the ratio remains at unity. In contrast to the case of only oscillations, this 
experiment is then unable to restrict m\. Not only that, since the argument [0] against the 
hierarchical four-neutrino spectrum depends crucially on the CDHSW experiment, it is clear 
that it cannot be valid in general. 

The idea of neutrino decay is of course not new. It is naturally related to the spontaneous 
breakdown of lepton number 0, |13|. The associated massless Nambu-Goldstone boson ||T4 



IS 



called the Majoron and the typical decay Majoron occurs if kinematically allowed. 

The triplet Majoron[0 is ruled out experimentally because the decay Z Majoron + 



partner (imaginary and real parts respectively of the lepton-number carrying scalar field) 
would have counted as the equivalent of two extra neutrino flavors. The singlet Majoron|^ 
is unconstrained because it has no gauge interactions. We assign lepton number L = — 1 
to z/s and assume the existence of a scalar particle x° with L = 2. [By convention, Vs is 
left-handed. If we use a right-handed singlet neutrino instead, then it would be assigned 
L = +1.] Hence the relevant terms of the interaction Lagrangian are given by 

Cint = gsl^sl^sX^ + hal^s{l^a(p° - L(P^) + h.C. (4) 

As and (0°) become nonzero, becomes massive and also mixes with i'e,fi,T to form the 
mass eigenstates 2^1,2,3,4- At the same time, \^Imx'^ becomes the massless Majoron M and 
the decay 

1^4 1^1,2,3 + M (5) 



is now possible. Neutrino decay involving only i'e,ti,T was recently proposed to explain 
the atmospheric data|l|], but that becomes a poor fit after the inclusion of the upward going 
muons[]TB[. More recently, it was shown [p!^ that combining oscillation and decay (at the 
expense of also adding Ug) gives again a good fit. In contrast, the effects we envisage here of 
z/4 decay in atmospheric and solar neutrino data are both small and do not change the usual 
oscillation interpretation appreciably, as shown below. 

Let i^e,fi,T,s be related to the mass eigenstates mi,2,3,4 through the unitary matrix Uai, 
which will be assumed real in the following for simplicity. Let » » m2 » mi 
with z/4 having the decay lifetime T4,. Then for solar and atmospheric neutrino oscillations 
with m\L/AE » 1, the probability of z/q, — is given by 

Pap = Sap{l - 2f/l) + Ul,U^p,{l + x^)-AY^ U^,Uo.jUp,Upj sin^ — (6) 

i<j<i 

where 

X = e-"'^^/^^^^ (7) 

In the case of laboratory experiments where AmfjL/4E << 1 for i < j < 4 but m\L/4:E is 
not necessarily large or small, the corresponding formula is 



0/3 



2Ul^ 1 - X cos ■ ^ 



2x cos — — + X 
2E 



2E 

Note that the above expression simplifies to a function of Uai, Up^, and x if is large, and 
to a function of UaA and Up4 alone if x = whatever the value of 1714. In those circumstances, 
the corresponding laboratory experiment has no sensitivity to oscillations, but does measure 
one fixed number. Specifically, if is large, then 

P,e = U'^4UU^ + x'), P^, = {I - Ul,f + x'Ut,, P,, = {l-Ul,f + x'Ul,. (9) 

If X = 0, then regardless of m4, Eq. (8) reduces to Eq. (9) but with x set equal to zero. The 
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LSND experiment obtains [I 



P^e = 3.1 ^1'^ ±0.5xl0-^ (10) 



whereas BNL-E734 has[|ig P^e < 1-7 x 10^^ and BNL-E776 haslTTJ P^e < 1-5 x 10 



-3 



Using the LSND 90% confidence-level limit of P^e > 1-3 x 10^'^, we find therefore reasonable 
consistency among these experiments. [The most recent result of the ongoing KARMEN II 
experiment ||18|| is P^e < 2.1 x 10"'^, which will eventually have the sensitivity to test Eq. (10).] 
The recent CCFR experiment [|12| measures P^e < 0.9 X 10~^ but its average L/E is one to 



two orders of magnitude smaller than those of the other experiments, hence its x-value may 
be taken to be close to one and the usual oscillation interpretation of the data holds. This 
constraint implies that m| < 30 eV^. 

At 1714 ~ 5 eV, we are below the CCFR exclusion and in a marginal region of the 
parameter space for pure neutrino oscillations consistent with the LSND evidence and the 
exclusion from BNL-E734 and BNL-E776. Between ~ 5 eV and ~ 3 eV, the BNL- 
E734 data exclude a solution if x = 1 and because that experiment has an average L/E an 
order of magnitude smaller than that of BNL-E776, LSND, or CDHSW, the decay factor 
goes against having a consistent solution here even if x < 1. Below 1714 ~ 3 eV, the oscillation 
+ decay interpretation of the latter 3 experiments becomes important, as shown below. 

Ideally, one should reanalyze the results of all the laboratory experiments using Eq. (8) 
and verify whether the positive LSND signal can coexist with the exclusion limits from the 
other laboratory experiments by extending the usual parameter space of m^, Ue^, and to 
include as well. This can be done only by using the full data set of each of the experiments 
and is best performed by the experimenters themselves. In the absence of such a calculation, 
we point out here the crucial fact that the CDHSW experimental would see no difference 
in its two detectors at distances of 130 m and 885 m, if the effective values of the quantity 
exp{—m4L/2ET4)cos{mlL/2E) is the same. In Table I, we show r4/m4(= l/T^mi) as a 



function of ml near 6 eV^ for which this happens, using as our very crude approximation 
the fixed values of Li/E = 0.065 m/MeV and L2/E = 0.442 m/MeV. This illustrates the 
possibility that the decrease from xi to X2 due to decay may be compensated by the increase 
in the value of the cosine from Li to L2 due to oscillations. Note also that there is a range 
of m| for which a null solution exists with varying T^/m^, whereas if the latter is zero, 
then m| has only discrete solutions (at 4.8 and 6.6 eV^ for example). In the realistic case 
of integrating over the experimental energy spectrum, both solutions will be smeared out, 
but the possibility of decay should result in a larger range of acceptable values of ml. For 
consistency, we also show in Table I the values of / = P^e/UlJJ'ji^ = l — 2x cos {m^L / 2 E) + x^ 
for the LSND and BNL-E776 experiments, using the fixed values of L/E = 0.75 and 0.5 
m/MeV respectively. This shows that the value of P^e as seen by the LSND experiment can 
be larger than that of BNL-E776 for 4.8 < mj < 5.8 eV^. 

To discuss solar and atmospheric neutrino oscillations, let us focus on the following 



specific model. Let cos^ = y2/3 and sin 6' = yl/S in Eqs. (1) and (2), and let mix with 
i>2 only, then Uai is given by 



U 



^2/3 







s^l/3 


-^1/6 




-\/l/2 


s^l/3 








s^l/3 





—s 





c 



where c and s are respectively the cosine and sine of the — 1/2 mixing angle. For solar 
neutrino oscillations, we have 

/ o2\2 . / 



COS 



Aml2L' 
2E 



+ 



x'^s'^ 
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(12) 



In the limit s = 0, this reduces to the usual two-neutrino formula with sin^2^ = 8/9 which 
is a good fit to the data0, either as the large-angle matter-enhanced solution or the vacuum 
oscillation solution. With a small of order a few percent [between 0.026 {x = 1) and 
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0.037 {x = 0) for P^e(LSND) = 1.35 x 10"^], this is definitely still allowed. Note that this 
result is not sensitive at all to the last term because s'^/Q is of order 10^^. 

For atmospheric neutrino oscillations, we have 

^ee=(l-y)' + ^, Pe. = P.e = (l+X^)^, (13) 

Here the limit s = corresponds to the canonical — > v^- solution with sin^ 29 — 1. 
As it is, the prediction of ~^ is still a fixed number, but smaller than unity (0.93 
for s^/3 = 0.037). Given that there is an uncertainty of about 20% in the absolute fiux 
normalization, we should consider instead the ratio 

where we have made an expansion in powers of and assumed that the ratio of to 
produced in the atmosphere is two. It is clear that this is numerically indistinguishable from 
the case s = . 

In this model, the decay z/4 — z/2 + ^ has some very interesting experimental con- 
sequences. For example, Ue from the sun decays through its component into i?2 = 
(c/ -\/3)(Pe + + ^t) ~ ^^s- Hence 

P{V, ^ Pe) = P{l^e ^ = P{^e ^ ^r) = ~ lO"', (16) 

where the energy of Pq is only 1/2 that of and x — Q has been assumed. This is in 
principle detectable especially since the v^.P capture cross section is about 100 times that 
of Vf.e scattering at a few MeV. Unfortunately, the Super-Kamiokande experiment has an 
energy threshold of 6.5 MeV for the recoil electron and taking into account the additional 1.8 

MeV threshold for the VeP e+n reaction, this would require the original Ve energy to be 
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2P -I- P 

Pee ~l~ '^P/j.e 



above 16.6 MeV, placing it outside the solar neutrino spectrum. With the recently lowered 
Super-Kamiokande energy threshold of 5.5 MeV, the fraction of solar Ve above 14.6 MeV 
is 1.6 X 10~^. Given the small probability of P(z^e ^e)i this will not change appreciably 
the total number of observed e-like events. Regardless of energy threshold, the inability of 
Super-Kamiokande to distinguish from e~ or to detect the 2.2 MeV photon from neutron 
capture on free protons makes it difficult to pin down this possibility in any case. 

In the Sudbury (SNO) neutrino experiment ||T9(] , the energy threshold for detecting recoil 



electrons is 5 MeV, but since there is also a threshold of about 4 MeV for breaking up the 
deuterium nucleus into two neutrons and a positron, the neutrino energy required is more 
than about 18 MeV. This again places it outside the solar neutrino spectrum. On the other 
hand, if the experimental energy threshold can be significantly lowered, then SNO may be 
able to see this effect because the z/g signature {i'e + d n + n + e^) is distinct from that of 

The best chance for detecting antineutrinos from the decay of z/4 is offered by the BOREX- 



INO experiment pO| with a very low energy threshold of 0.25 MeV. Taking into account the 
1.8 MeV needed for inverse beta decay, i.e. v^P ~^ e~^n, this means that solar neutrinos with 
energy above 4.1 MeV can be detected as antineutrinos. The idea of looking for antineu- 
trinos from the sun was motivated by the possibility of a large neutrino magnetic moment 
which may convert Ue into z/g in the sun's magnetic field. The capability of BOREXINO for 
detecting this has been discussed earlier ||2l|. For our new distinctive effect of 1/4 decay, the 



observed antineutrino energy spectrum is predicted to go from f{E) to f{E/2), where E is 
the energy of the original neutrino. 

For atmospheric neutrinos, since and z/g are produced together with z/^ and z/g in about 
equal amounts, it is not possible to tell if a given event comes from the primary neutrino or 
its decay product, even if the detector could measure the charge of the observed lepton. 
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To search for the Ue transition in the LSND and KARMEN experiments, one 

would use the monoenergetic (29.8 MeV) z/^ from 7r+ decay at rest, which has the signature 
of a monoenergetic positron of 13.1 MeV from inverse beta decay, i.e. UeP — > e~^n, in 
coincidence with a 2.2 MeV photon from the subsequent capture of the neutron by a free 
proton. However, this signal is overwhelmed by the neutral-current reaction u ^^C —>■ v ^^C*, 
with the subsequent emission of a 15.1 MeV photon. 

In proposed long-baseline z/^ appearance experiments, the oscillation probability is 
given by 

which is not easily distinguished from the s = case. However, the decay products of 1/4, 
i.e. z/g, t'^, and z/^, may be observable with their own unique signatures, depending on the 
capabilities of the proposed detectors. 

In the case of four-neutrino oscillations, the effective number of neutrinos in Big Bang 
Nucleosynthesis is an important constraint p2[ . In this model, with m4 ~ few eV and ~ 



few percent, the presence of a stable would have counted as an extra neutrino species, 
making = 4. This may not be acceptable if < 4, as indicated from the observed 
primordial ^He abundance |2^. The decay of z/4 changes to 3 + the contribution of the 



Majoron (i.e. 4/7). 

With z/4 as a component of z/g, neutrinoless double decay has an effective mass of 
(s^/3)m4 ~ 0.2 eV if ~ 5 eV. This value is just at the edge of the most recent experimental 



upper boundp4|. 



Finally a comment on the neutrino contribution to dark matter may be in order. With z/4 
decaying and mi, m2, and being too small, there is no neutrino dark matter. However, 
it is possible that mi ^ m2 ~ ~ few eV, while m4 is higher by another few eV, in 
which case ui, z/2, and z/3 will contribute to dark matter. Our discussion goes through almost 
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unchanged, except that ml in Eq. (8) will be replaced by m| — 2,3- 

In conclusion, we have shown in this paper that a hierarchical four-neutrino scenario 
is acceptable as a solution to all present neutrino data regarding the disappearance and 
appearance of Ve and f^. The assumed singlet neutrino of a few eV may decay into a 
linear combination of the three known doublet neutrinos with half of the energy. This 
new feature allows our proposal to be tested in future solar neutrino experiments such as 
BOREXINO (and perhaps SNO), and should be considered in forthcoming long-basehne 
accelerator neutrino experiments. 
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/(LSND) 


/(E776) 


1.8 





3.92 


0.01 


5.0 


0.030 


3.01 


0.01 


5.2 


0.065 


2.21 


0.19 


5.4 


0.085 


1.72 


0.38 


5.6 


0.095 


1.37 


0.57 


5.8 


0.095 


1.09 


0.78 


6.0 


0.086 


0.54 


1.03 


6.2 


0.068 


0.55 


1.37 


6.4 


0.038 


0.22 


1.94 


6.6 





0.0 


3.0 



Table 1: Null solution for oscillation and decay at the two CDHSW detector distances. 
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